gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十四章/eg14_13.m
wc2=1.5; num=20; den=conv([1 1 0],[1 2]); G=tf(num,den); [mag,phase,wcg,wcp]=margin(G); margin(G) t1=1/(0.1*wcg); beta=10; Gc_lag=tf([t1,1],[beta*t1,1]) G1=G*Gc_lag; [mag,phase,w]=bode(G1); mag1=spline(w,mag,wc2); L=20*log10(mag1); alfa=10^(-L/10); t2=1/wc2/sqrt(alfa); Gc_lead=tf([alfa*t2,1],[t2,1]); G0=G*Gc_lead*Gc_lag; figure(2) margin(G0) figure(3) step(feedback(G0,1))